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Abstract 

Renormalization factors for AS = 1 four-quark operators in the effective 
weak Hamiltonian are perturbatively evaluated in domain-wall QCD. The 
one-loop corrections of AS" = 1 four-quark operators consist of two types of 
diagrams: one is gluon exchange between quark lines and the other is penguin 
diagrams containing quark loops. Combining both contributions, our results 
allow a lattice calculation of the amplitude for K — ► tttt decays with 0(g 2 ) 
corrections included. 
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I. INTRODUCTION 



AS = 1 four-quark operators appear in the effective low-energy Hamiltonians for non- 
leptonic weak decays of kaons, which are relevant for the AI = 1/2 rule and the CP- 
violation parameters e'/ e. Although comparison of their experimental results with theoretical 
predictions based on the standard model is supposed to give a good testing ground for 
the model, it is not accomplished without any reliable non-perturbative estimates for the 
hadronic matrix elements of AS = 1 four-quark operators. 

Lattice QCD calculation allow us to evaluate the AS = 1 hadronic matrix elements 
from first principles in QCD. While in the past decade the Wilson and the Kogut-Susskind 
quark actions are used to calculate the four-quark hadronic matrix elements, these quark 
actions have inherent defects: the explicit chiral symmetry breaking in the Wilson quark 
action, which causes the non-trivial operator mixing among different chiralities, and the 
mixture of space-time and flavor symmetries in the Kogut-Susskind quark action, which 
annoys us by demanding non-trivial matching between continuum and lattice operators. The 
domain-wall quark model [jlJQ is a five-dimensional Wilson fermion with free boundaries 
imposed on the fifth direction of N s layers. This quark formulation is expected to have 
superior features in the large N s limit: no fine tuning for the chiral limit, 0(a 2 ) scaling 
violation and no mixing between four-quark operators with different chiralities. Perturbative 
calculations have shown that these expectations are fulfilled at the one-loop level in the 
limit N s — > oo PH71 , and numerical studies support these features non-perturbatively 



These advantageous features over other quark formulations urge us to apply the domain-wall 
quark to a calculation of the AS" = 1 hadronic matrix elements. 

Since the Wilson coefficients in the effective low-energy Hamiltonians for non-leptonic 
weak decays of kaons are calculated in some continuum renormalization scheme (e. g., MS), 
the corresponding AS* = 1 hadronic matrix elements should be given in the same continuum 
scheme, which requires us to convert the matrix elements obtained by lattice simulations to 
those defined in continuum renormalization scheme. A necessary ingredient for this trans- 
formation is the renormalization factors matching the lattice operators to the continuum 
ones. 

In this article we make a perturbative calculation of the renormalization factors for 
the AS" = 1 four-quark operators consisting of physical quark fields in the domain-wall 
QCD(DWQCD). The one-loop corrections of AS = 1 four-quark operators contain the 
gluon exchange diagram and the penguin diagram, which leads to the following form for the 
relation between the lattice and continuum operators: 

QT nt = E ztiQf + z f cn Qltn + o( 9 % (l) 

3 

where Qi denote a set of AS" = 1 four-quark operators. The gluon exchange contributions Zfj 
are already calculated in our previous paper , whose results are applicable to AS* = 2 four- 
quark operators. In this work we evaluate the additional contribution denoted by Zf en Q 1 ^ 
which originates from the penguin diagram in the AS = 1 case. With the aid of the results 
in Ref. 0, we obtain complete expressions of the renormalization factors for the AS* = 1 
four-quark operators. 

This paper is organized as follows. In Sec.|T| we present the calculations of the penguin 
diagram in the continuum and on the lattice. Full expressions of the renormalization factors 
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for the AS" = 1 operators are given in Sec.|T|. Our conclusions are summarized in SecpV] 
In Appendix we give the DWQCD action and the Feynman rule relevant for the present 
calculation to make this paper self-contained. 

The physical quantities are expressed in lattice units and the lattice spacing a is sup- 
pressed unless necessary We take SU(iV) gauge group with the gauge coupling g, while 
N = 3 is specified in the numerical calculations. Hereafter, repeated indices are to be 
summed over unless otherwise indicated. 



II. PENGUIN DIAGRAMS 
A. AS—1 four-quark operators 

We consider the following AS = 1 four-quark operators, 



Qi = 


(s a U b ) V - A {Ubda)v-A, 




(2) 


Q 2 = 


{s a Ua)v-A{u b d b )v-A, 




(3) 


Qs = 


(s a d a ) V -A ^2(qbqb)v-A, 

Q 




(4) 


Q 4 = 


{s a d b ) v _ A ^2(q b q a )v-A, 
q 




(5) 


Q 5 = 


{s a da)v-A y^X<lbqb)v+A, 
1 




(6) 


Q 6 = 


{s a d b )v-A y^Xqbq a )v+A, 
q 




(7) 


Qr = 


3 

- (S a d a ) V-A J2 e<? V- 

1 q 


rA, 


(8) 


Qs = 


3 

-(s a d b ) V -A e9 fe<?a)v- 

1 q 




(9) 


Q9 = 


3 

-{s a d a ) v _ A ^e q {q b q b ) v . 
1 q 


-A, 


(10) 


Qio = 


3 

-0a4)v-A e<1 (Ma)v- 

1 q 


-A, 


(11) 



where V ± A indicates the chiral structures 7^(1^=75). a, b denote color indices and e q are 
quark charges: 2/3 for the up-like quarks and —1/3 for the down-like quarks. The flavor 
sum is intended over those which are active below the renormalization scale of the operators. 
This set of operator basis closes under QCD and QED renormalization. We classify these 
10 operators into the following four types: 



Qi = {s a r x u b ){u b r x d a ), (12) 

Q2 = (s a T x u a )(u b T x d b ), (13) 

Q\ = (s a T x d b ) J2 <4(q b T£q a ) i = 4, 6, 8, 10, (14) 

q 

QI = (s a r x d a ) «< (ftlyft) 2 = 3,5, 7, 9, (15) 
q 
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where = 7 M (1 - 75) and T Y = 7 M (1±7 5 ). aj are given by 



a\ = \ i = 3,4,5,6, (16) 
ot q i=\e q 1 = 7,8,9,10. (17) 



B. Continuum calculation 

It is instructive to first show the calculation of the penguin diagram in the continuum 
regularization schemes. For the present calculation we employ the Naive Dimensional Reg- 
ularization(NDR) and the Dimensional Reduction(DRED) as the ultraviolet regularization 
scheme, in each of which the loop momenta of the Feynman integrals are defined in D 
dimensions parameterized by D = 4 — 2e (e > 0). The major difference between both reg- 
ularization schemes are found in the definitions of the Dirac matrices: NDR scheme defines 
the Dirac matrices in D dimension, while DRED scheme retains them in four dimensions. 

We calculate the following Green function with massless quarks: 



{Qi)ab;cd = (QiSa(Pl)d b (P2)q C (P3)qd(P4)) , (18) 

where a, b, c, d label color indices, while spinor indices are suppressed. Truncating the exter- 
nal quark propagators from (Qi), where we multiply ip'j (j = 1,2,3,4) on (Qi), we obtain 
the vertex functions, which is written in the following form up to the one-loop level 



Mp)} 



ab:cd 



Af(p)+Ar sj (p)+Ar pj (p) 



(1-9) 



(1-P), 



ab:cd 



(19) 



where the momentum p is defined by p = p\ — p 2 = P4 — p%. The superscript (i) refers to 
the i-th loop level. At the one-loop level Aj- 1-9 ^ and Aj 1 p ' represent the gluon exchange 
diagram and the penguin diagram, respectively. In this paper, however, we do not consider 
A^ 1 ^ because the gluon exchange contributions to the renormalization factors are already 
calculated in our previous paper |J. 

The tree-level vertex functions A 4 - ^ are given by 



Aj 0) 

e(0) 
i 

o(0) 



ab;cd 
ab;cd 
ab;cd 
ab\cd 



101 
l®i 



ab.cd 



1 X 



(X- 



ab.cd 

101 



1 x 



p/Li 



p/Li 
J- \ 



101 



ab\cd 



ab\cd 



X 



X 



X ■ 



1 Yi 



p/Lt 

1 Yt 



(20) 
(21) 
(22) 
(23) 



where represent the tensor structures in the Dirac spinor space and © and 
color space as 



act on the 



101 

'mi 



ab:cd 



= Sad5, 



ad^cbi 



ab:cd 



= SnhS, 



abOcd- 



(24) 
(25) 
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We take q = u in eq.fllq) for convenience. This special choice, however, does not affect the 
final results for the renormalization factors. 

The penguin diagram is illustrated in Fig. [I], which yields the one-loop level vertex 
functions A,- 1- ^, 



A (i-P) 
At P) 



A 



e(l-p) 



A 



o(l-p) 



ab:cd 



ab:cd 



ab:cd 



ab:cd 



0. 



-9 



+9' 



T A ®T B 
T A ®T B 



T ®T 



B 



G AB ( P ) L aP (p) tr(T s Ylal ^) I* x ® 7v , 



(26) 
(27) 
(28) 

lu, (29) 



where T A (A = 1, . . . , A^ 2 — 1) are generators of color SU(AQ. G AB is the gauge propagator 
given by 



G A J(P) 




NDR, 
DRED, 



(30) 



with D-dimensional metric tensor 5^ u and the four- dimensional one 5^. L a p{p) denote the 
integral of the quark loop, whose result is 



L a (p) 



6(4tt) 
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P 



- + \og(fi 2 /p 2 ) + - + 1 ^5 af3 +[- + \og{n l /p 1 ) + - PaPp 



(31) 



where 1/e = 1/e — 7 + ln |47r| and /i is the renormalization scale. We remark that the metric 
tensor 5 a p is defined in D dimension both for the NDR and DRED schemes. The products 
of the Dirac matrices are reduced as follows: 



tr [T Y jaltilp) <W = 4 x 
tr [Ty7a7/*7/9] PaPf3 = 4 x 

TxlctliflfFyPcVfi = S X y x 



(2 - D)5, 



fJtU 

28^ - DS^ 



NDR, 
DRED, 

f 2 PtlPv - b^p 2 NDR, 
1 2 PtM p u - ~5^p 2 DRED, 

' 2(L> - 2) 2 T/i (l - 7s ) NDR, 
-4 [27^ - D lfl ] (1 - 75) DRED, 

2{D - 2) [p 2 lfl - 2pJ] (1 - 75) NDR, 
4 [p 2 lfl - 2pJ] (1 - 75 ) DRED, 



(32) 
(33) 
(34) 
(35) 



where 7^ = in the DRED scheme, and we choose = 7^ — -^7^ as an evanescent 

operator in the DRED scheme [11]. It should be noted that the terms proportional to p 1 are 
eliminated by imposing on-shell conditions on external quark states. 

After some algebra we obtain the following expressions for the vertex functions (p6|) — (p9|) : 



At P) 



ab;cd 



12tt 2 



T A ®T 



ab:cd 



7m(! -75) ®7m L 



cont 



(36) 



5 



A e(l-p) 
t 

Ao(l-p) 



ab:cd 



ab;cd 



12tt 



E 



a' 



T ®T / 



ab\cd 



TmC 1 -7s) ®7 m 



12tt 2 



T A ®T A 



ab\cd 



7m(1 -7s) ®7m ^ c 



where 



with 



I + log(/i 2 /p 2 ) -h I + Cj 



Q 2 and Q° in NDR, 
Qf in NDR, 

Q 2 , <?? and Q° in DRED. 



(37) 
(38) 

(39) 



(40) 



The pole term 1/e is subtracted from Lf ont in the MS scheme. 

From the tree-level vertex functions in eqs. (|20"D — (|23|) we find 



T A ®T' 



ab:cd 



7 (U ( 1 -7s) ®7 M 



101- —1(8)1 
Af + A< 0) - - 



7m(! - 75) ® - 75) + TmC 1 + 75)] 

ab;cd ^ 



N 



(A 



(o) 
3 



A 



(oh 



Using this relation we finally obtain 
(1-p) _ 9 2 1 



A 



12vr 2 4 



C(Q0 L 



cont 



Ai 0) + Ar-^(Ai 0) 



(41) 



(42) 



with C(Q 2 )=1, C(Q 3 ) = 2, C(Q 4 ) = C(Q 6 ) = f q , C(Q 8 ) = f u - f d /2, C(Q 9 ) = -1, 
C(Qio) = fu — fd/2 and C(Qi) = for other z, where / g , and denote the number of 
flavors, the number of charge 2/3 up-like quarks and the number of charge —1/3 down-like 
quarks, respectively. 

C. Lattice calculation 

Let us turn to the calculation of the vertex functions on the lattice. In this subsection we 
use the same notations for quantities defined on the lattice as those for their counterparts 
in the continuum. 

We consider the Green function of eq. flilf) on the lattice. At the one-loop level the 
penguin diagram in Fig. [T] contributes to the vertex functions: 

A (i- P ) 



A 



(i-p) 



A 



e(l-p) 



A 



o(l-p) 



ab:cd 



ab:cd 



ab:cd 



ab:cd 



0. 



-9 



+9' 



-9 



T A ®T B 
T A ®T B 



T ®T 



B 





(1- 


. ab;cd 




(1- 


ab;cd 






(1- 


- ab;cd 



1 " <?G A *{p) £ aftr(T« L»)ri ® 7,, 



Gfiip) [afrlL„(p)r£ + a°T 5 Y L,{p)T 5 x 



(43) 
(44) 
(45) 

(46) 
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where (1 — Wq) 3 is the overlap factor to the four-dimensional quark fields which emerges 
through the truncation of the external quark legs by multiplying zj^ (i = 1, ... ,4) on the 
Green function. denotes the integral of the quark loop, which is expressed by 



LM = £ —{P R S F (k) ls + P L S F (k) NsS }(-i) Vti (k - p/2) st 

x {S F (k - p) tl P L + S F (k - p) tNs P R } (47) 

(48) 



with 



I^k,p) = -c fl (k-p/2) i _^ {a+a , ) lfl - c^k - p/2) FFI ^ _ e _ {a+a>)) farf 



+ir Sfl {k - P/ 2 ) F / (1 _ e -(a + aO) ^ + iTS ^ k ~ P/ ^J{Y 



-(a+a 1 ) 



(49) 



where c^(k) = cosfc M , s At (fc) = sinA;^, s'^ = s^{k — p) a = a(k), a' = aik — p), F = F(k) and 
F' = F(k — p). The expressions of v^, S F , a and F are given in Appendix. 

Since the function L^(p) has an infrared divergence for p 2 — > 0, we consider extracting 
its divergent part by employing an analytically integrable expression I^(k,p) which has the 
same infrared behavior as /^(/c,^), 

r* d 4 k ~ „ x r d A k 



LM = £ 7^yi4(M + £ 7^)4 - 4(M} , (50) 

where we choose 



fcM^A'-^l-^,^, (51) 

with A (< 7r) a cut-off. Now the second term in the right hand side of eq.(|50|) is regular in 
terms of p, we can make a Taylor expansion, 

r* d 4 k ~ „ , r* d 4 £; 



+ L j^l pppx dp^p- x ^ p) ~ ~^ p) } (52) 

where p, A = 1,2,3,4. Before evaluating each term of the right hand side, let us consider 
the transformation properties of L M under various discrete symmetries, from which we can 
predict the types of terms to be allowed in L M . The same kind of discussion is found in the 
Wilson quark case [|T2|]. 

Under the operation of charge conjugation matrix C = 7472 we find 
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Cv^{k) st C 1 = -vT(-k) Ns+1 - t>Ns+1 - s , (53) 

CSpikUC- 1 = Sp(—k) N a +i-t,N 3 +i-s, (54) 
CL,{p)C~ l = -Lj(p). (55) 

Parity transformation gives 

7 4 f 4 (A;4, fc)st74 = ^(^4, -k)N s +\-s,N a +i-t, (56) 

^4Vi(k 4 , k) s tlA = -Vi(k 4 , -k) Ns+ i_ StNs+1 - t , (57) 

j 4 S F (k 4 , ^)si74 = S F (k 4 , -k) Na+1 ^ SjNs+1 ^ t , (58) 

74^4(^4,^)74 = L 4 (p 4 , -p), (59) 

74^i(P4,p)74 = ~Li(p 4 , -p), (60) 

where i — 1, 2, 3. These discrete symmetries restrict the form of L M as 

with cx M! , = py M , 7^/2, where the massless quark case is considered. We can eliminate further 
terms by the Ward-Takahashi identity: 

2sin(p M a/2)^(fc - p/2) st = [S F {k)-% t - [S F (k - p)" 1 ]^, (62) 

which leads to 

2sin(p M a/2)L M (p) = i J^-0^[P R {S F (k) u - S F (k - p) n } P L 

+P R {S F (k) 1Ns - S F (k - p) 1Ns }P R 
+P L {S F (k) Nsl - S F (k - p) Nsl } P L 
+P L {S F (k) NsNs - S F (k - p) Ns n s }Pr] 
= (63) 

where the periodicity of S F is used. Under the requirement that the terms of the left hand 
side must vanish order by order in terms of the lattice spacing, we find Cq = c 2c = and 
c 2a = —C2b in eq.(|6l|), with which the expression of is simplified as 

Lfi{p) = -cxia^pv + c 2 a{l^p 2 - Pvft) + 0{a). (64) 

(X 



Here we should note that 



7<5<V7<5 = 0, (65) 
tr ((7^(1 ±75)) = 0, (66) 



which means a^ u p v term in eq.(|64D does not contribute to the wave functions (|43D — (f46|) . 

Now we see that L^O) = and linear term in p is irrelevant in the expansion (|52"D, we 
focus on the remaining terms. Performing the integration of I^k^p) — I^k, 0) analytically, 
we obtain 
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r d 4 k 



{i,(k, P )- 1,(^0)} =(i- 



1 



tL'r 



J-n (2tt) 4 I ^ ^' ^ ' 7 J v u '16tt 2 

As for the last term in the expansion of eq. (|5^) 
expression for the integrand: 

d 2 



(~P\+Pj) I (log(A 2 /p 2 ) + 5/6) - jl» 
some tedious algebra leads to the following 



P P P\ 



dppdpx 



{l - I n(p)} \ P =o 



~^~Pfi1fi (.Rfi Rfifi ~ Rfijtv "I" — + Rfxfjtfj, + R^vu ~ Rfifi^v ~ Rfi^vft 

+lvl a PvPalniRi^ v ~ Krfv^a), 



where R£, R^, R^ and B* VQ are given by 



r; 



R 



— e" 



4(e c 



+ 



s 2 f 



c M (10s 



4F 2 (1 -e 



-2a N 



5S M 



4F(e° - e" 

,2 



„ -e-) 2 
■V 15 - 1 a* fjy* fin 



rs. 



+ 



2rs 4 .e a 



4F(e Q - e" a ) F(e 
32(1 - Wq)0(A 2 - k 2 



a _ e -a\3 



flie a + e~ a ) - g^(e a - e 



I.-' 



(e Q - e" Q ) : 

Cfi4 



(A; 2 ) 4 ' 

s 2 J 2 Ae a + e~ a ) - (h u + s 2 v g uu )(e a - e~ a ] 
(2s 2 -s 2 )(H u + slG uu ) 



F 2 (l - e 



-2a 



+ C, 



+re- a s 2 JslG uu + H v ) 



(68) 



(69) 



R 



flV 



R 



'fiua 



+re c 



?-a\3 



F(e 

-{I - wl)6{X 2 - k 2 ) 



[slo™ + K)(e a 



slflie a + e 



, „ t , ;2N 2(P-2A; 2 )(4fc 2 -^ 



(k 2 y 



2F 2 (t 



2 9 2 2fi 
I 2 _/£ I fi fi v fi v 
+ S v?v p + F 



2F{e a 



+ re 



S fl S U^IJ,v S U C ^R^ ^^yRu/ 1 ^ 



+re c 



C fi S vfv 



2F(e Q 



-a\2 



+ re c 



s 2 s 2 



F(e a 

+ (1 - w 2 )6(A 2 - k 2 )l6k 2 uJ Sy v 



-a\3 



9fivi< 



(70) 



-2 



C ^ F " Ca -4(l- W 2 )fl(A 2 -A^ ^ 



(P) 3 ' 



(71) 
(72) 
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with 



— rW + c M + r cosh a 

W sinh a 
—c u rW + cf. — s 2 , + rc u cosh a 



W sinh a 



9l 



a' 



9u^ 



ffifu 



+ 



r 



tanh a W sinh a 
2a -r + Wf, f, 
11 C F 2 (e a -e- a ) 



_l_ r 



w ' 



F„ = e 



i-' 



c 2a rc^-Wh^ | 



F(e a - e~ a ) 2 ' 
h„ 



F 2 (e a - e~ a ) F(e a - e~ a ) 



(73) 
(74) 
(75) 
(76) 

(77) 



G 



G Vil - F 



9\xu 



-a\2 



e a + e 



+ e' 



M-r + WfJ + M-r + Wf, 



-a\2 



F 3 (e a - e' 



2e 2a ( 



F 2 (e a - e 

r + WfJi-r + Wf v ) + e a F{-r(U + f») + WfJ u + Wg, v } 



(75 



We choose r = — 1 in this calculation. Using the results of eqs. (|67|) and flBlf ), L^p) is 
expressed as 



:i-^o 2 ) 

167T 2 



p\) \ (log(A 2 /p 2 ) + 5/6 



where 






A = 


1 167T 2 


1 r 


6 1 — vol 


2J- 


B = 


16vr 2 1 r* 


d 4 k 




(27T)' 


AR = 


167r 2 1 r 71 




1 - Itfg 2 i-TT 


(2tt)« 


5R = 


167T 2 1 




1 - Wl 2 J-n 


(2tt)< 



(79) 



/9 7? c _|_ f? 1 ^ _|_ I?^ I 
( jd(1 rid \ — n 

4 y-^nau -^nvaj u > 



(80) 
(81) 
(82) 
(83) 



with jj, ^ v 7^ a. We find that 5R = from the symmetry of the integrand. The integrals 
are estimated by a mode sum for a periodic box of a size L 4 with L = 64 after transforming 
the momentum variable through p^ = — sing^. We choose A = 7r for the cut-off in the 
integrand 1^. The numerical results for A are B are presented in Table p] as a function of 
M. We observe that the expected relation A = — B is well satisfied. We also checked that 
AR is consistent with zero as expected. Finally we obtain 



'l-wp L latt 
16tt 2 3 



(84) 



where 
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L latt = log(7r 2 /p 2 ) + ~ + SB. 





Substituting the above expression for L^(p) in eqs.([4l3l) — (f46|), we have 



Ai 



(i-p) 





ab;cd 




pY 




9 2 




ab;cd 


127T 2 


pY 




9 2 




ab\cd 




pY 


ab;cd 


= 5 X Y . 



2\4 r latt 



1 - 



</ /' lal I 

2\4/„ s 



ab;cd 

T A ®T A 



7^(1 -75) ®7 M , 

7m(!-75) ®7**j 



c(\ r latt 



ab;cd 

T A ®T A 



(85) 

(86) 
(87) 
(88) 



ab:cd 



7m(1-75)®7m> (89) 



where we used the on-shell conditions for the external quark states to eliminate the terms 
proportional to In the same way to derive eq.([4"2"D in the continuum calculation, the lattice 
wave functions are written in the compact form as follows, 



A 



(i-p) 



9 2 1 



12tt 2 4 



C{Qi) L 



latt 



Af+A^-^Af +Af) 



(90) 



where C(Qi) are already obtained in the previous subsection. 

Comparing the continuum vertex functions in eq. ([42]) and lattice ones in eq. (|9~0"D, we find 
that the penguin diagram contributions to the renormalization factor in eq.(|l]) are written 

as 



1 



1 



zr n Q l ;ii = ^w zr [Qa +q "- n {Q3 + g5)]latt ' 

where the penguin operator and its coefficient are given by 



(91) 



Q 



latt 
pen 



latt 



^pen 



g 2 C(Qj 

12vr 2 4 



latt 



L 



cont \ 



g 2 C(Qj 

16tt 2 3 



— log(/ia) 2 + z, 



pen 



(92) 
(93) 



with zf = log(7r ) + 3B q. Numerical values of zf with the DRED scheme are 

6 

presented in Table | function of M. 



III. FULL RENORMALIZATION FACTORS 

In order to write down the complete expressions for the renormalization factors of the 
AS = 1 operators Qi, we still need to know the contributions from the gluon exchange 
diagrams, which is denoted by i^Q' att in eq. (|l|). Fortunately, they can be obtained from 
the combinations of the results in our previous paper |J, where we calculated the gluon 
exchange diagrams for the following four-quark operators, 
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0±(qi, 92, 93, 94) = ^ [{.qiq2)v-A{.qzqA)v-A ± {qm)v-A{q3q2)v-A 
1 



^1(91,92,93,94) 
^2(91,92,93,94) 



-CF(9i92)y-A(9394)y+A + (giT g 2 )v-A(93^ 94)v+a 
1 



2iV 



(9i92)y-A(9394)y+A + (qiT g 2 )y-A(93^ 94)y+A 



(94) 
(95) 
(96) 



with Cf = (N 2 — l)/(2iV) the second Casimir of SU(iV) group. The tensor structures in the 
color space are 1®1 and T A ®T A . 

The operators Qi are related to C±,i,2 as 



Q 2 ±Qi 


= 80±(s, u, u, d), 


(97) 




= 8'£0±(s,d,q,q), 
q 


(98) 


Q9 i Qio 


= 8]T— (9 ± (M,9,9), 
q 1 


(99) 




= 8j20 2 (s,d,q,q), 

q 


(100) 


Qs 


= 8 ]T— 2 (M,9,9), 
g 1 


(101) 


-NQ 5 + Q 6 


= 8j20 1 (s,d,q,q), 

q 


(102) 


-NQ 7 + Q s 


3p 9 

= 8£^-0i(M,?,<z)- 

q 1 


(103) 



From these relations the renormalized operators, which include the contributions of both 
the gluon exchange and penguin diagrams, are expressed in terms of the lattice operators 
as follows: 



Q 



cont 



(i - <?zi 



i^g £-jlatt _|_ ^penQlati 



where 



Zl = { 



1 + 



1 + 



1 + 



9* 



16tt 2 
_l_ 
16tt 2 
_l_ 
16tt 2 



^log(/ia) 2 + ^p 
3 

-— log(/ia) 2 + zi- v 2 i 

3(iV 2 -l), . , 2 

— log(//a) + z 2 + v 2 i 



pen 



i = 1,2,3,4,9, 10, 
i = 5,7, 
% = 6,8, 



(104) 



(105) 



7 9 



7 9 



16tt 2 



-31og(/m) + 



2 , z + ~~ 2 - 



16tt 2 



3 log(yua) + 



2 , ^2 - zi + u 2 i - «12 



TV 



16tt 2 



Nv 



21 



(i,j) = (l,2),(3,4),(9,10), 

= (5,6),(7,8), 

(2,j) = (6,5),(8,7), 
others 



(106) 
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with £±,1,2) V12, and V21 presented in our previous paper ||. The penguin operator Q 1 ^ 
and its coefficient Zf en are given in eqs. (^) and fl93|) , respectively. For the sake of conve- 
nience, we list the differences between the NDR and DRED schemes in the finite part of the 
renormalization constant, 



y zi - v 2 i 



y z 2 + v 2 i 



z 2 - z 1 + v 2 i - v 12 \ 



1 NDR 



1 NDR 



NDR 



NDR 



N 



-) 

/^pen\NDR 



Z+ + Z_ 



(zi - V21) 
{z 2 + V 2 1 



DRED 



DRED 



.DRED 



N 2 -6 



2N 



N 2 



2N ' 
N 2 -4 



DRED 



5 
2' 



z 2 - Zj + ^21 - f 12 
iV 

iV( U21 ) DRED - 3, 
5 

(2 pen )D RED + ' 



DRED 



7 
2' 



- z 



- z 



2, 3, 5, 7, 9, 



4,6,8,10, 



(107) 

(108) 

(109) 

(110) 

(111) 
(112) 

(113) 



where vi 2 = v 21 = in the DRED scheme. 

We are interested in the magnitude of the renormalization factors with the currently 
accessible coupling constant in numerical simulations. Let us estimate it taking f3 = 6.0 
with M = 1.8 as a representative case. All the necessary ingredients in this analysis are 
given in our previous papers |5]|| and Table |T] of this report. With the use of the mean field 
improvement |13], we have 



(l-w (M) 2 )Z™ F (M) 
z™ F (M) + z™ F (M) 



V MF 



(M) 



,MF 



(M) 



zf F {M) = -11.187, 
z f F (M) = -18.6, 
zf Y {M) - zf F {M) -- 



= 0.9085, 
-12.848, 

2.491, 
-7.4, 



.pen 



(M) = 2.328, 



in the DRED scheme, where we employ 



9m0-/a) 
M 



P- -0.13486 
6 

M + 4(n- 1) = 



= 2.1792, 
1.3112, 



(114) 
(115) 

(116) 

(117) 
(118) 
(119) 
(120) 



(121) 
(122) 



with P 



u 



0.59374 the plaquette value at (3 = 6.0 in the quenched approximation. 



Combining these results we obtain 
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MF 



u 2 x 0.9968 1 = 1,2,3,4,9,10, 

u 2 x 1.0245 1 = 5,7, (123) 
u 2 x 0.9006 z = 6,8, 

z g xMF r = Z% = 0.0416 (z,j) = (l,2),(3,4),(9,10), 

-0.0412 (z,j) = (5,6),(7,8), (124) 

others 



:i-^o 2 ) 2 ^. 

^pen 

:i-i) 2 ^. 



MF 



C(Qi) x 0.0130, (125) 



at fia = 1 in the DRED scheme, where we factor out the tadpole contributions. We find 
that the penguin diagram contributions to the renormalization factor are quite small. 

IV. CONCLUSION 

In this paper we have presented the full expressions for the renormalization factors of 
the AS = 1 four-quark operators up to the one-loop level including the contributions of 
both the gluon exchange and penguin diagrams. Our calculation, however, does not include 
mixing with lower dimensional operators sd and s^^d. The coefficients of these operators 
diverge as inverse powers of the lattice spacing, due to which it is practically inappropriate 
to subtract these lower dimensional operators by the perturbation theory. We are instead 
forced to use the non-perturbative methods such as those based on the chiral perturbation 
theory jl4| and the Schrodinger functional scheme Hl5"| . We leave it to future work. 
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APPENDIX 

In this appendix we explain the domain-wall fermion action and its Feynman rules rele- 
vant for the present calculation. The domain- wall fermion action is written as [0, 



N s 

Sow = X 



n s =l 



2 



+ - U{n) s {l + 7 5 )^(n) s +i + ^(n).(l - 7s)V>W s -i) + (M - 1 + 4r)^(n) s ^(n), 



2 

+ mJ2 {^Wn.PrWti)! +^{n) l P L iP{n) Ns ) , (126) 

n 

where is the link variable of the SU(iV) gauge group and the Wilson parameter is set 
to r = — 1. Four dimensional space-time coordinate is labeled by n and s is an extra fifth 
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dimensional index which runs from 1 to N s . Since we impose no gauge interaction along the 
fifth dimension, it is also possible to consider that the index s labels the N s "flavor" space. 
In our one-loop calculation we take N s — > oo to avoid complications arising from the finite 
N s such as mixing among the four-quark operators with different chiralities. The parameter 
m denotes the physical quark mass and at m = one chiral zero mode is supposed to appear 
under the condition < M < 2 for the Dirac "mass" M. Pr/l are projection operators 
defined by Pr/l = (li 7s)/2. For the gauge part we employ a standard four dimensional 
Wilson plaquette action. 

The quark fields on the four-dimensional space-time are given by the combinations of 
the fermion fields at the boundaries, 

q{n) = Ptfl>(n)i + P^p(n) N „ (127) 
q(n) = ^{n) N P R + $(<n) x P L . (128) 

We consider the composite operators constructed with these "physical" quark fields, and our 
renormalization procedure is based on the Green functions consisting of only these quark 
fields. 

In order to obtain the Feynman rules we perform the weak coupling expansion of the 
quark and gauge actions. The fermion propagator with momentum p is obtained by inverting 
the domain- wall Dirac operator in eq . f p. 26| ) , which is expressed by S F (p) s t as an N s x N s 
matrix in the "flavor" space. In the present one-loop calculation, however, we do not need 
the whole matrix elements because we consider the Green functions consisting of the physical 
quark fields. The relevant fermion propagators are restricted to following three types: 

(q(-p)q(p)) = PrS f (p) u Pl + P L S F (p) NsNs P R + P L S F (p) Nsl P L + P R S F (p) 1Ns P R 
—27^ sinp^ + (1 — We~ a ) m 
= - (1 - e a W) + m 2 (l - We~ a ) ' 

?(-p)^«(p)) = PrSf(p)u + PlS f {p) Ns s 

^ (i^sm Ptl - m (l - We- a )) (V^^Pr + e~ a ^P L 



F 

1 

+ F 



with 



+- (e- a{s ~ 1] P L + e- a{Ns - s) P R ) e~ a [m (i 7 „ sinp^ - m (l - We 



W = l-M-r$^(l-cosp M ), (132) 

n 

cosh(a) = — 7| ' ( 133 ) 

F = 1 - e a W - m 2 (l - We~ a ) , (134) 



(129) 



m (i7 M sin^ - m (l - We' 01 )) - P]e~ Q (e- a(s ~ l) P R + e- a{Ns - s) P L ) , (130) 
(ip s (-p)q(p)} = S F (p) sl P L + S F (p) sNs P R 
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where the argument p in the factors a and W is suppressed. 

In the perturbative calculation of Green functions we assume that the external quark 
momenta and masses are much smaller than the lattice cut-off. In this case the external 
quark propagators can be expanded in terms of them. We have the following expressions as 
the leading term of the expansion: 



l-wl 



{qq)(p) = - 

if + (1 — WQ)m 

= (qq)(p) (<~ 1 Pl + Wo 

(M(p) = (w^Pr + wZ'-Pl) (qq)(p), 



,N s ~sp 



(135) 

(136) 
(137) 



where wo = 1 — M. The form of (qq) (p) tells us that y 1 — Wq is the overlap factor to the 
four-dimensional quark fields at the tree-level. 

The one-gluon fermion vertex with outgoing and incoming fermion momenta p and q, 
respectively, is given by 



V 1 1(q,p)st = gT A v,( 



q^ + Pv 



St 



gT A {r^ cos (^^) + r sin (H^) R " (138) 



where T A (A = 1, . . . , N 2 — 1) are generators of color SU(iV). We do not need two-gluon 
fermion vertex in the present calculation. 

The gluon propagator for a gluon of momentum p is written as 



Gi B ( P ) 



>AB 



4sin 2 (p/2) 



<V - (i - « 



4sin(p M /2) sin(p„/2) 
4sin 2 G»/2) 



(139) 



where sin 2 (p/2) = sin 2 (p M /2) and we choose the Feynman gauge (a = 1) in our calcula- 
tion. 
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TABLES 



TABLE I. Numerical values for A, B and z? en (DRED) as a function of M. Note that zf en is 
independent of % (% = 1, . . . , 10) in the DRED scheme. 



M 


A 


B 


z pcn (DRED) 


0.05 


-1.5078(4) 


1.5077(3) 


5.729(1) 


0.10 


-0.1130(4) 


0.1128(3) 


1.544(1) 


0.15 


0.2836(4) 


-0.2838(3) 


0.355(1) 


0.20 


0.4456(4) 


-0.4458(3) 


-0.131(1) 


0.25 


0.5193(4) 


-0.5195(3) 


-0.352(1) 


0.30 


0.5513(4) 


-0.5515(3) 


-0.448(1) 


0.35 


0.5607(4) 


-0.5609(3) 


-0.476(1) 


0.40 


5565(4) 

V / • '.7 ' J \J ' J \ X / 


-0 5567(3) 


-0.464(1) 


0.45 


0.5434(4) 


-0.5436(3) 


-0.425(1) 


0.50 


0.5242(4) 


-0.5244(3) 


-0.367(1) 


0.55 


0.5003(4) 


-0.5006(3) 


-0.296(1) 


0.60 


0.4728(4) 


-0.4729(3) 


-0.213(1) 


0.65 


0.4419(4) 


-0.4421(3) 


-0.120(1) 


0.70 


0.4081(4) 


-0.4082(3) 


-0.019(1) 


0.75 


0.3712(4) 


-0.3714(3) 


0.092(1) 


0.80 


0.3314(4) 


-0.3316(4) 


0.211(1) 


0.85 


0.2884(4) 


-0.2885(3) 


0.341(1) 


0.90 


0.2417(4) 


-0.2419(3) 


0.480(1) 


0.95 


0.1911(4) 


-0.1913(3) 


0.632(1) 


1.00 


0.1360(4) 


-0.1362(4) 


0.797(1) 


1.05 


0.0760(11) 


-0.0754(4) 


0.980(1) 


1.10 


0.00879(7) 


-0.00944(65) 


1.178(2) 


1.15 


-0.0648(3) 


0.0647(4) 


1.400(1) 


1.20 


-0.1476(3) 


0.1472(4) 


1.648(1) 


1.25 


-0.2407(4) 


0.2405(3) 


1.928(1) 


1.30 


-0.3470(4) 


0.3468(4) 


2.246(1) 


1.35 


-0 4695(4) 


4693(4) 


2.614(1) 


1.40 


-0 6125(4) 


6122(3) 


3.043(1) 


1.45 


-0 7818f4) 


7816(4) 


3.551(1) 


1.50 


-0.9859(4) 


0.9857(3) 


4.163(1) 


1.55 


-1.2369(4) 


1.2367(3) 


4.916(1) 


1.60 


-1.5531(4) 


1.5529(3) 


5.865(1) 


1.65 


-1.9638(4) 


1.9636(3) 


7.097(1) 


1.70 


-2.5179(4) 


2.5177(3) 


8.759(1) 


1.75 


-3.3049(4) 


3.3047(4) 


11.120(1) 


1.80 


-4.5053(4) 


4.5051(3) 


14.721(1) 


1.85 


-6.5451(3) 


6.5449(3) 


20.841(1) 


1.90 


-10.7173(1) 


10.717(3) 


33.357(1) 


1.95 


-23.549(7) 


23.549(7) 


71.853(22) 
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FIGURES 



Pi P 2 




FIG. 1. The penguin diagram, s, d and q label quark flavors, pi, p2, P3, Pa, P and k are 
momenta. Solid circle denotes the four-quark operator. 
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